Structural anomalies at the magnetic transition in centrosymmetric BiMnC>3 
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The structural properties of BiMn(93 were determined from neutron powder diffraction data as a 
function of temperature and magnetic field. The structure at all temperatures was found to be cen- 
trosymmetric with space group C2/c, which is incompatible with ferroelectricity. At Tc ~ lOOif , we 
observed the onset of a large magnetoelastic strain, proportional to the square of the magnetization. 
We interpret this structural rearrangement, together with the previously observed magnetodielectric 
anomalies, as due to the need to optimize the partially frustrated magnetic interactions. 

PACS numbers: 75.25,+z, 75.30.Kz, 75.80.+q, 77.80.-e, 77.80.Bh, 
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Multiferroic materials, for which two or all three of 
ferroelectricity, (anti-) ferromagnetism and ferroelastic- 
ity are observed in the same phase, have received re- 
newal interest in recent years. Such systems are rare 
in nature but are potentially interesting for a wide ar- 
ray of technological applications 0, 0]. In particular, 
BiMn03 has been extensively studied in both bulk and 
thin film forms as a candidate ferromagnetic (FM) fer- 
roelectric (FE) material. BiMnOs shows FM ordering at 
the Curie temperature Tc = 99-105 K, a unique feature 
among the Y*>\BOz transition- metal series (B — Cr-Ni). 
The observed ferromagnetism is related to ordering of 
the partial filled 3d 2 2 and 3& x 2_ y 2 orbitals on Mn, as 
evidenced by a cooperative Jahn- Teller distortion taking 
place well above Tc [3|]. Dos Santos et al. showed that 
BiMn03 is monoclinic above and below Tc, and pro- 
posed the non-centrosymmetry space group C2, which in 
principle allows ferroelectricity along the b axis [4( . How- 
ever, establishing experimentally that BiMn03 is also FE 
has proved considerably more difficult, due to the lack of 
single crystals and the low resistivity of BiMnOs films. 
Band structure calculations performed in C2 by Hill and 
coworkers [5j predicted the existence of ferroelectricity in 
BiMnOs, and Shishidou et al. later calculated a polar- 
ization of 0.52 C/cm 2 [6|. dos Santos et al. [3] showed 
FE hysteresis loops at room temperature and below, on 
both bulk and thin films samples, with a FE polariza- 
tion of 0.043 C/cm 2 at 200 K, one order of magnitude 
lower than predicted. Also, Sharan et al. reported ev- 
idence for second-harmonic generation (SHG), strongly 
enhanced by applied electric field, consistent with the 
presence of FE domains in their thin- film samples [8]. All 
these data would indicate that BiMnOs is a proper FE, 
since ferroelectricity is observed already at room temper- 
ature. This rules out any mechanism based on symmetry 
breaking at the onset of magnetic ordering, as recently 



proposed, for example, for TbMnOa[9]. In support of 
these results, a separate set of measurements shows quite 
convincingly a very strong magneto-electric coupling in 
BiMnOs near Tc- Kimura et al. [10[ investigated the 
magnetodielectric effect in BiMnOs bulk samples: they 
observed a pronounced anomaly in the dielectric constant 
e, and a very strong magnetic field dependence of e near 
Tc, which they ascribe to the coupling between FE and 
FM order. 

This substantial body of work has been very recently 
challenged by Belik et a/. [HI] After a detailed analy- 
sis of convergent beam electron diffraction and neutron 
diffraction data on bulk samples, Belik concluded that 
there is no evidence for inversion symmetry breaking in 
BiMnOs, and proposed the centrosymmetric space group 
C2/c, which cannot support ferroelectricity. It is note- 
worthy that first-principle calculations on BiMnOs also 
show that the centrosymmetric C2/c structure is more 
stable than the C2 structure at K [l2[. The simplest 
explanation for this clear contradiction between these 
data and previous literature has to do with sample dif- 
ferences, particularly between films and bulk. In par- 
ticular, the lower Tc of thin films with respect to the 
bulk may indicate partial oxidation, leading to different 
polymorphs [13[ . Thin films are also affected by the pres- 
ence of the substrate, which could reduce their symme- 
try. However, Kimura's data [10] are impossible to dis- 
miss on this ground, since they were collected on high- 
quality polycrystalline sample. Nevertheless, the obser- 
vation of a large magnetodielectric effect does not require 
ipso facto the breaking of inversion symmetry, and may 
arise from structural anomalies of a different kind. Here, 
we present neutron diffraction measurements collected on 
BiMnOs as a function of temperature (10 < T < 295 K) 
and magnetic field (0 < H < 10 T). In agreement with 
Belik et a/ [ill]. we find that the centrosymmetric space 
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FIG. 1: Main panel: Magnetization vs. temperature for 
BiMnOs, collected on both heating and cooling in a 100 Oe 
magnetic field. The anomaly at 42 K (arrow) is due to a 
small amount of antiferromagnetic Bi2Mn40io, which does 
not show hysteretic behavior. Inset: Mn magnetic moment 
versus temperature, as determined from neutron diffraction 
data (see text). 



group C2/c is perfectly adequate to model our data at 
all temperatures and fields. We observed very significant 
structural anomalies upon entering the magnetically or- 
dered phase. On the basis of the Landau theory of phase 
transitions, we propose that the large magnetodielectric 
effect observed by Kimura is related to the anomalies 
we observe, most likely through a subtle change in the 
environment of the highly polarizable Bi ion, and both 
are due to the need to optimize the partially frustrated 
magnetic interactions. 

Polycrystalline samples of BiMnOs were synthesized 
by solid state reaction of Bi2 03 (Aldrich 99.99 %) and 
Mn 2 3 (Aldrich 99.999 %) at 40 Kbar and 600 °C for 
3 hours, using a high-pressure multi anvil apparatus. 
Three samples were synthesized at the same experimental 
conditions in order to obtain a suitable quantity of pow- 
der for neutron diffraction experiment (~ 2g). Magne- 
tization curves were measured in the temperature range 
5-20 K in steps of 5 K and from 20 to 110K in steps 
of IK with a superconducting quantum device (SQUID) 
magnetometer (Quantum Design). The magnetization 
curves, measured in a 100 Oe magnetic field both on 
cooling and in heating mode (Fig. [1]), evidenced the pres- 
ence of a small thermal hysteresis (Tc (heating) = 100.26 
K, Tc (cooling) = 99.05), consistent with the magnetic 
transition having a slight first-order character. 

Neutron powder diffraction data were collected at 
room temperature (RT) using the high resolution powder 
diffractometer (HRPD) and as a function of temperature 



using the general materials diffractometer (GEM) at the 
ISIS facility. Data were collected on warming from 10 
to 293 K in steps ranging from 1 to 10 K using a he- 
lium cryostat. Different collection times (20 min- 6 h) 
were employed to obtain higher statistics at key tem- 
peratures. Additional measurements were performed on 
GEM at a fixed temperature (95K) using a cryomagnet 
to vary the magnetic field from to 10 T, in steps of 0.5 
T. Rietveld refinements of the structure (nuclear above 
Tc and nuclear and magnetic below Tc) were performed 
at all temperatures in both C2 and C2/c space groups, 
using the program GSAS In the course of the re- 

finements, the magnetic moments on the inequivalent Mn 
sites were constrained to be equal, and to lie along the 
high-symmetry y direction. 

The structural parameters as refined on the RT HRPD 
data in C2 and C2/c are reported in Tab. I of the Sup- 
plementary Material. Both the quality of the refinements 
and the actual structural parameters are extremely sim- 
ilar in the two cases, resulting in similar bond distances 
and comparable formal oxidation states for Bi and Mn. 
In spite of the larger number of parameters, the R factors 
obtained in C2 are slightly higher than those in C2/c. In 
addition, the estimated standard deviations of the atomic 
coordinates are one order of magnitude larger in C2, indi- 
cating the presence of significant correlations. We there- 
fore conclude that the centro symmetric space group C2/c 
is the appropriate one to describe the structure. The 
same conclusion applies to the structural and magnetic 
refinements of the 10K data collected on GEM. A com- 
parison of the RT and 10 K structures (Tab. II and III 
of the Supplementary Material) indicates that the dis- 
tribution of shorter-basal and longer-apical bonds is es- 
sentially the same, confirming that the orbital ordering 
pattern is unchanged on cooling. The refined magnetic 
moment per Mn site as a function of temperature is re- 
ported in the inset of Fig. [IJ showing a good agreement 
with bulk magnetization measurements. 

The pattern of orbital ordering giving rise to ferromag- 
netism is best visualized by referring it to the pseudo- 
cubic perovskite axes (Fig. [2]), which also make it very 
easy to compare BiMnOs with other orbitally-ordered 
perovskites such as LaMn03. The appropriate cell is 
actually triclinic, it encompasses 4x4x4 primitive per- 
ovskite units and is related to the monoclinic cell by the 
covariant transformation 



[ftp, bp, Cp] [&m? t) m , C m ] 



with 



4 x 



A, 



' 1 


1 


1 


2 





-2 


_ 1 


-1 


1 


b p 




4 



(1) 



3.98 A, 



3.92 

ol p = 7p = 91.4° and j3 p = 90.8°. Fig. [2] shows a section 
of this structure perpendicular to the c p axis. Subsequent 
layers are staggered by one perovskite unit along the a p 
direction. One should also bear in mind that a n — b n 



3 




FIG. 2: color online. Orbital ordering pattern in BiMnOs. 
Main Panel: A slice cut perpendicular to the pseudo-cubic 
c p axis, showing Mn06 octahedra with in-plane (red) and 
out-of-plane (light blue) long Mn-0 bonds, corresponding to 
occupied d z 2_ r 2 e g orbitals. Inset: the same structure shown 
on the monoclinic axes. The presentation emphasized the 
occupied orbitals. 



and c p — bp sections are identical, since they are related 
by the mirror symmetry in C2/c. In Fig. O dark red 
octahedra have their partially filled orbital (long bonds) 
within the projection plane, whereas light blue octahe- 
dra have it perpendicular to it. In the case of LaMnOa, 
the same section (through the b axis in the Pnma set- 
ting) would only show dark red octahedra with identi- 
cal arrangement. Based on the well-known Goodenough- 
Kanamori rules, it is therefore easy to understand that 
the dark red "stripes" carry a uniform FM interaction 
as in the equivalent LaMnOs layers. The magnetic in- 
teractions along the b p direction are also uniformly FM. 
Interactions around "light blue" octahedra are partially 
frustrated, 4 FM and 2 antiferromagnetic (AFM), but the 
FM interactions, which are also expected to be larger in 
magnitude, clearly dominate, giving rise to an overall FM 
structure. 

The bottom and middle panels of Fig. [3] display 
the evolution of the monoclinic lattice parameters as a 
function of temperature. The dimensionless normaliza- 
tion ((p(T) - p(RT))/p(RT)) for p = a, b and c and 
/3(T) — /3(RT) for the monoclinic angle in radiants) is 
appropriate for the calculation of the strain tensor (see 
below). There is a clear anomaly with onset at Tc for 
the monoclinic angle and, to a lesser extent, for a, while 
b and c are hardly affected by the transition. The solid 
lines represent fits to the data above the transition us- 
ing a simple Einstein model of the phonon thermal ex- 
pansion, with T#=200 K (the results are affected only 
slightly by the exact value of T#). This enabled us to 
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FIG. 3: Magneto-elastic effect in BiMnOs- Bottom: Rel- 
ative change ( J^sook) ) °^ ^ ne m onoclinic lattice parameters 
p = a, b and c and Middle: change of the monoclinic (3 an- 
gle vs. temperature. The parameters were determined from 
Rietveld refinements based on neutron diffraction data (sym- 
bols). The lines represent the calculated thermal expansion 
due to phonons, based on a simple Einstein fit to the high- 
temperature data with Te — 200 K(see text). Top, excess 
thermal evolution (i.e., with the phonon contribution sub- 
tracted off) of the pseudo-cubic lattice parameters. Lines are 
guides to the eye. If not shown, the error bars are smaller 
than the symbols. 



extract the " anomalous" component of the lattice distor- 
tion due to magnetoelasticity (see below). The top panel 
of Fig. [3]shows this anomalous distortion after projection 
onto the pseudo-cubic axes (the effect of the transition 
on the triclinic angles is very small). This presentation 
clarifies the origin of the magnetoelastic effect: the lat- 
tice shortens along the unfrustrated direction b p (Fig. 2, 
inset), and it lengthens to a greater extent along the two 
other partially frustrated directions. The overall effect on 
the unit cell volume is a slight expansion on cooling (neg- 
ative magnetostriction, not shown). The magnetoelastic 
origin of these lattice anomalies finds further support in 
the data we collected as a function of applied magnetic 
field at fixed temperature (T = 95K, just below Tc). Fig. 
[H demonstrates that the anomalous variation of the (3 
angle (with the small phonon component subtracted) is 
proportional to the square of the ordered magnetic mo- 
ment, regardless of whether the moment is induced by 
temperature or field. 

In order to relate these changes to the magnetodielec- 
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FIG. 4: Change of the monoclinic f3 angle as a function of 
the square of the ordered magnetic moment on the Mn site, 
as obtained by varying the temperature in zero field (circles) 
and by varying the applied magnetic field in the range 0-10 
T at a constant temperature of 95 K (squares). 

trie anomalies observed by Kimura et al. [10[ , we need to 
consider an appropriate form of the Landau free energy 
including the lattice strain. Since the structure remains 
monoclinic below Tc, the strain does not break the sym- 
metry, and all the elements of the strain tensor are in 
principle allowed to contribute independently to the free 
energy. However, for simplicity, we will only consider the 
scalar strain e sc , essentially the quadrature sum of the 
elements of the strain tensor. In the monoclinic case, it 
can be shown [l5[ that only the 633 and e3i tensor ele- 
ments contain A/3, and that, for the relevant values of the 
parameters and ignoring the small a anomaly, e sc ~ nf3, 
with n ~ 1. Therefore, Fig. [H provides a quantitative es- 
timate of the magnetoelastic strain. A similar result can 
be obtained by considering the diagonal elements of the 
strain in the triclinic setting (Fig. 02 top). The magne- 
toelastic strain is one order of magniture larger than the 
one reported by Kimura on the basis of dilatometric mea- 
surements - a discrepancy most likely due to the powder 
averaging in the latter. In the Landau free energy expan- 
sion, the terms responsible for the magnetodielectric and 
magnetoelastic coupling are |M 2 P 2 and ^e sc M 2 + f e 2 sc , 
respectively, where M is the magnetization, P is the po- 
larization and and c are phenomenological constants. 
By minimizing the free energy, one can show that, near 
the transition and for large values of the dielectric con- 
stant e, 

e sc = --M 2 (2) 
c 

In other words, both the dielectric anomaly and the 



magnetoelastic strain are proportional to the square of 
the magnetization, and therefore to each other. However, 
since the coupling constants involved are different, there 
is no reason why the proportionality constant should be 
of the order 1 (it is in fact ~ 20). Here, we have clearly 
shown the presence of a structural rearrangement upon 
magnetic ordering in BiMn03, as evidenced by very sig- 
nificant lattice anomaly. However, the most likely micro- 
scopic origin of the observed dielectric anomalies is not 
the change in lattice metrics itself, but rather a related 
subtle change in the environment of the highly polariz- 
able Bi 3+ ion. Our data both above and below the mag- 
netic transition confirm that, within our accuracy, the 
structure of BiMn03 retains a center of inversion, and is 
therefore very unlikely to be ferroelectric. 
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